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Abstract 

We show that large- A" <j) A theory is not trivial if one accepts the presence of a tachyon with a truly 
huge mass, and that it allows exact calculation. We use it to illustrate how to calculate the exact 
resummed pressure at finite temperature and verify that it is infrared and ultraviolet finite even in 
the zero-mass case. In 3 dimensions a residual effect of the resummed infrared divergences is that at 
low temperature or strong coupling the leading term in the interaction pressure becomes independent 
of the coupling and is 4/5 of the free- field pressure. In 4 dimensions the pressure is well defined 
provided that the temperature is below the tachyon mass. We examine how rapidly this expansion 
converges and use our analysis to suggest how one might reorganise perturbation theory to improve 
the calculation of the pressure for the QCD plasma. 



1 Introduction 

In traditional calculations of the pressure in massless quantum field theories at finite temperature, 
there is a breakdown of perturbation theory because of infrared problems [1] [2]. This is true, in 
particular, of gauge theories at finite temperature. Formally the same types of problem crop up also 
in scalar \<j) A theory. The traditional procedure [1] for summing the ring diagrams (figure la) produces 
a series expansion for the pressure P in powers (and logarithms) of A 1 / 2 . This takes care of infinitely 
many otherwise problematic higher-loop diagrams, but it does not cover diagrams that include a 
single one-particle-irreducible self-energy insertion and so leaves arbitrarily high-loop orders to be 
considered. In nonabelian gauge theories, where some propagators remain massless perturbatively, 
the two-particle- irreducible higher-loop diagrams are still potentially infrared divergent. In an earlier 
paper [3] we have proposed a different and less direct method which absorbs all of the diagrams into 
a single resummed one-loop quantity which enjoys manifest infrared regularity in four dimensions. 
This introduces a variable mass m for the field as a parameter, and gives the pressure in terms of an 
integral over the thermal propagator corresponding to mass m. 

In this paper, we apply our summation method to sum all the extended ring diagrams of figure lb, 
which we call foam diagrams. Other authors [4] [5] have called these super-daisy or Hartree-Fock 
diagrams. We began this analysis as a warm-up exercise for the gauge-theory problem, but have found 
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Figure 1: (a) ring diagrams; (b) foam diagrams 
that it has interesting features in its own right. 

Moreover, it allows us to investigate the convergence properties of thermal perturbation theory in a 
solvable model, from which we derive suggestions on how to optimize the perturbative treatment of 
the pressure also in more complicated theories such as QCD. 

If we consider an O(N) scalar field theory with unrenormalised Lagrangian density 

C(x) = \ ({d$(x)f - mMx)f) - ((fe) 2 ) 2 (1.1) 

then in the limit N — > oo the pressure per scalar particle coincides with the one obtained from the 
infinite sum of foam diagrams in the N = 1 theory. Hence our calculation in this paper of the foam 
diagrams may be regarded either as the leading-TV term in the pressure in the large-iV theory, or 
as an approximation (we do not know how good) to the pressure in the N = 1 theory. We obtain 
an expression for this infinite sum that is derived from a path integral and so has validity beyond 
perturbation theory. The renormalisation of the mass, and in 4 dimensions also of the coupling, is an 
essential part of the analysis; this is achieved through Dyson equations which can be solved exactly, 
and so again goes beyond perturbation theory. 

We describe the renormalisation in section 2. As is well known [6], and as we find, in 4 dimensions 
there are problems with 4 theory. For reasonably small coupling, these problems turn out to arise 
only at huge mass scales and so are not really important for physics, but in order to handle them 
we write our equations for n dimensions. In section 3 we review our formula [3] for the pressure; as 
before, we choose to use the real-time formalism of thermal field theory. The formula involves the 
thermal addition 5m 2 to the renormalised squared mass m 2 . A key simplification is that, for the foam 
diagrams, 5m 2 is both independent of momentum and real, and we show how this leads to rather 
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simpler versions of our general formula for the pressure. We need to solve an integral equation for 
(5m 2 , which we do in section 4 for the case of 3 dimensions. We show that the infrared divergences 
that occur in the zero-mass case are rendered harmless by our resummation and that in the strong- 
coupling or low-temperature limit the leading term in the interaction pressure becomes independent 
of the coupling and is 4/5 of the free-field pressure. 

In section 5 we investigate what happens in 4 dimensions if we ignore the triviality difficulties of the 
theory. We find that the integral equation for 5m 2 has no solution above some critical temperature 
T ma ,xi so that it seems that then the pressure does not exist. However, for sufficiently weak coupling, 
T max is exponentially huge, so that the 4-dimensional theory may be accepted as a highly accurate 
effective theory. In the particularly interesting massless case we make a detailed comparison of our 
nonperturbative results with those of resummed perturbation theory. We find that the latter converges 
quickly only for rather small coupling and that the rate of convergence depends critically on the 
renormalization scale. We end with some speculation on the case of QCD, where explicit calculation of 
the pressure up to and including order g 5 has revealed rather bad convergence properties of resummed 
perturbation theory there. We suggest that it may well be fruitful to reorganise the expansion so that 
Sm/T, rather than the coupling g, is the expansion parameter. 



2 Renormalisation 



Zero temperature 

We choose a zero-temperature renormalisation scheme that makes the formula for the pressure P(T) 
as simple as possible. In lowest-order perturbation theory the renormalised mass is defined to be 

m 2 = m 2 + \ M(ml) (2.1a) 

where M corresponds to the single-loop Feynman graph: 



M(m 2 ) = \\ -^-- - (2.2a) 

V ' 2 I {2ir) n q 2 - m 2 + it 



We extend this to all orders of perturbation theory by replacing (2.1a) with the Dyson equation 

m 2 = m 2 + \ M(m 2 ) (2.16) 

The integration over q° in the integral in (2.2a) may be done by closing the contour in one or other 
half plane and taking the residue at the pole. The result is that instead we may write 

M(m2) = 3 / J^ 27T6+iq2 ~ ^ ^ 

The integration is simple: 

M (m 2 ) = V 2 J m n - 2 (2.2c) 

V ' 2(4vr)™/2 v > 

When we go to 4 dimensions we need also to renormalise the coupling. Initially we keep n / 4 as a 
regulator. We define the renormalised coupling A to be the value of the %i — > jj scattering amplitude 
at s = 0, where i and j denote "colour" labels. The equation for A is shown diagrammatically in figure 
2. If we work to leading order in N in the large- ./V theory, or choose to include only foam diagrams in 
the N = 1 case, we must omit the last two terms in figure 2, so that we have 

A = A + A AL(m 2 ) 
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Figure 2: coupling renormalisation 

L(m 2 ) = \\ S^TT-o % ^ = M'(m 2 ) (2.3a) 

y ' 2 J (27r) n (g 2 - m 2 + «e) 2 v ' y ' 

where the prime denotes differentiation with respect to m 2 . Both M and M' are ultraviolet divergent 
when n — > 4. The coupling has dimension m 4_n . From (2.3a) and (2.2c), 

An 



1 + C„A m"- 4 
T(2 - ^n) 

C n = \ - V 2.4a 
2(4yr) n /2 ^ ; 

so that if we want both the bare and the renormalised coupling to be non-negative 

< A < — - (2.5) 

Hence, when n — > 4, A vanishes for all values of m, because C n diverges. This is the well-known 
triviality of c/> 4 theory in 4 dimensions [6]. 

One might perhaps believe that it does not matter whether the bare coupling is positive, because 
only the renormalised theory is relevant. However, if for n = 4 we choose A to be greater than 0, the 
renormalised theory has a tachyon. To see this, we write the equation corresponding to figure 2 (with 
the last two terms omitted because they are negligible in the large- iV limit) for the i% — > jj scattering 
amplitude for a general value of s: 

T(s) = X + X M'(m 2 ,s)T(s) 



M'(m 2 ,s) 



From (2.4a) and (2.6), 



1 f d n q i 

2 J (2tt)™ ((q + \p) 2 -m?+ ie) ((q - \p) 2 - m 2 + it) 

-C n [ dx (m 2 - sx(l - x)) n/2 ~ 2 (2.6) 
Jo 



T(s) = j (2.7a) 

1 + C n X Jo dx{(m 2 - sx(l - x)) n/2 ~ 2 - m"" 4 } 



which becomes in 4 dimensions when \s\ is large 



1 - 32^ log(-s/m 2 ) 
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It is evident that, in 4 dimensions, T(s) has a pole at 



2 /32tt 2 \ , , 

s = stachyon ~ -m exp I — ^— I (2.8a) 

Our definition (2.3a) of the renormalised coupling A in terms of the value of the scattering amplitude 
T(s) at s = makes A vary with the renormalised mass m, and the appearance of a tachyon pole in 
T(s) is related to the occurrence of a Landau pole in A(m 2 ). Although this tachyon pole is in principle 
unacceptable, in practice it occurs at a very large negative value of s. For example, if A = 1, it is at 
— 10 137 m 2 . One could avoid the presence of this tachyon by introducing an UV cut-off A, since then 

A = (2.46) 

1+ 3^1og(A 2 /m 2 ) 

which makes it possible to have both A and A positive, with A < 327r 2 / log(A 2 /m 2 ) and therefore 

I ^tachyon | ^ A . 

Or alternatively we simply say that, for reasonably small coupling, the tachyon is so far away that it 
can be ignored. Either way, the n = 4 theory seems to be perfectly acceptable as an effective theory, 
for energies much below either the tachyon mass or the cut-off. 

In what follows we shall be particularly interested in the case of massless theories, because there one 
encounters the infrared problems of thermal perturbation theory. Because of the tachyon (2.8a), this 
case may seem to be excluded, for then the tachyon cannot be kept far away, unless we send A — > as 
m — > 0. Contrary to appearances, this does not signal a trivial theory, however. In the massless case, 
there are infrared divergences in the defining equation for the renormalised coupling (2.3a). Switching 
to an alternative scheme defined by 

A = A + A AL(/2 2 ) (2.36) 
with some arbitrary scale p, shows that then, when A / 0, 

T(s) = = (2.7c) 

l-3A,[log(- S /A 2 )-2] 

vanishes at s = 0, but not otherwise. So whenever we are going to inspect the massless limit, we shall 
switch over to using A rather than A. In terms of A and p, the tachyon is at 

■Stachyon = ~p 2 exp + (2-86) 

so for given choice of p we can always make it as far away as we wish by taking A small enough. 

In minimally subtracted massless theories, a renormalization scale is usually introduced through A — > 
fi 4 ~ n X. Because 

,4-n r /^2\ 1/2 i__ P- 



^- n L(^)~-— - _- 7 -lo; 



32vr 2 \A-n 1 & 4vr^ 2 

the scale p, 2 introduced as in (2.3b) is seen to coincide with the one of the modified minimal subtraction 
scheme (MS), p 2 = 47re~ 7 ^ 2 . Hence our notation. 

A more physical renormalization scheme would-be to define the renormalized coupling through the 
scattering amplitude at some scale Jl 2 so that A = T{—]1 2 ). (2.7c) shows that p 2 = e 2 p 2 . However, 
in conformity with most of the literature on thermal field theory, we shall continue to use the MS 
scheme. 



5 



Nonzero temperature 

At nonzero temperature, in addition to the renormalised zero-temperature mass there is a thermal 
contribution, because (2.1) becomes replaced with 

m 2 + 5m 2 =ml + X M T (m 2 + 5m 2 ) 

M T (m 2 ) = M{m 2 ) + N T {m 2 ) 
2, f d n q , 2 2 1 



NT{m ) = J (2# 2 ^ 2 " m2) ^7T3T ( 2 - 9 ) 
Eliminating the unrenormalised mass m\ with (2.1b) gives 

5m 2 (m 2 ,T) = X [M T (m 2 + 5m 2 ) - M(m 2 )] (2.10) 

For less than 4 dimensions both Ao and the expression in square brackets are finite, but when n — > 4 
this is no longer true and we must introduce also the renormalised coupling. If we define the function 

M(m 2 ,5m 2 ) = M(m 2 + 5m 2 ) - M(m 2 ) - 5m 2 M\m 2 ) 

= C n {{m 2 +5m 2 ) n/2 - 1 - (m 2 ) n ' 2 - 2 {m 2 + (n/2 - l)5m 2 )} (2.11a) 

we find that 

5m 2 (m 2 ,T) = X[M(m 2 ,5m 2 ) + N T (m 2 + 5m 2 )] (2.12) 
This is an integral equation for 5m 2 (m 2 ,T) in which everything is finite when n — > 4, with 

M(m 2 , 5m 2 ) -^r{{m 2 + 5m 2 ) log ( 1 + - 5m 2 \ (2.116) 

32ir 2 I \ m l ) J 

The expansion of this in powers of 5m 2 jm 2 begins with a term proportional to (5m 2 ) 2 . In perturbation 
theory 5m 2 = 0(A), so that AM = 0(A 3 ) and the first term on the right-hand side of (2.12) represents 
a contribution that first appears at three loop order. Since frequently calculations up to two loops 
have been taken as a clue to exact results, this particular contribution has been repeatedly missed 
in the literature [4] [7]. It owes its existence to the interplay between zero-temperature contributions 
plus their renormalisation with the thermal effects, and it will turn out to be of great importance to 
the existence and behaviour of the solutions of the equation for 5m 2 . 

In the massless limit m — > 0, we are again confronted with infrared divergences, which can be avoided 
by switching to the alternative renormalisation scheme (2.3b) This amounts to substituting 

which, in the limit m — > 0, leads to 

5m 2 {0,T) = X[M{fL 2 ,5m 2 ) + N T (5m 2 )] (2.14) 

with 

I 5m 

M(fl 2 ,5m 2 ) = ^Sm 2 (log - 1) (2.15) 

It is easy to see that 5m 2 (0,T) is, as it should be, independent of p, once the /x-dependence of A is 
taken into account. 
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3 The pressure 



Our formula [3] for the pressure is derived from the grand partition function 

Z(T) = j d$ exp{iS[faT}) 

S[faT}) = [ d n xC(x) (3.1) 
Jc 

In the variant of the real-time thermal field theory we use* , the integration is over all x and over the 
time contour C which, in the complex t-plane, runs along the real axis from — oo to +oo, back to — oo, 
and then down to — oo — i/T. (This is known as the Keldysh contour [9].) Differentiate with respect 
to ttiq keeping Ao fixed: 

d 



dm 2 



log Z = -\iZ~ l J d$ (^j d n x$ 2 {x)^j exp(iS[faT}) 

= -\i{ [ d n x^{x)) (3.2) 
Jc 

Here, < . . . > denotes a thermal average. Space-time translation invariance tells us that the thermal 
average of 4> 2 {x) is independent of x, and so the x integration is trivial: 

° log^ = -^(<? 2 (0)> 



dm 2 6 2T 



NV f d n q 19/ . . 

]J2^ D ^ T) (3 - 3) 



2T 

Here, 



D (q, T) = d n x e^ x fa (0)fa (x) (3.4) 



where fa denotes any component of the 0(iV)-symmetric field <p. D 12 (q, T) is an element of the familiar 
2x2 matrix propagator D(q, T) of Keldysh-contour real-time thermal field theory [2] [8]. Inserting its 
known form [10] [11] into (3.3) and using the relation P(T) = (T/V) logZ, we obtain 

dP ( T ) - f ^ Im 1 f35a) 

dm 2 J (2tt)" e« a /T _ i f _ m 2 _ n ( g , T, m) K j 

When, as in the application to the foam diagrams, the self energy n is real, we must apply the usual 
it prescription to m\ and so 



9P(T ) _ ,.. 2 
dm 2 



M T (m 2 + U(q,T,m)) 

o 

2 , r2\ 



= M T (m 2 + 5m 2 ) (3.56) 

where Mt is defined in (2.9). Subtracting off the zero-temperature pressure and using (2.10), we have 

d /„ /rTnN „, \ \ 5m 2 (m 2 ,T) 
M (P(T)-P(0))= L_lJ (3.6) 

* For a review of the essentials of real-time thermal field theory, see reference [8] 
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Using (2.1a) and (2.3a) to express everything in terms of renormalised quantities and imposing the 
obvious boundary condition that the pressure vanishes for infinite mass, we find the remarkable formula 

, . r° , /2 5m 2 (m' 2 ,T) . „ , 

p(T) - F(0) =i/ m (3 ' 7a) 

Here, 

A(m ' 2) = l + A^-^K 2 )] {m 
so that A(m' 2 ) is a running coupling that is equal to A when m' 2 = m 2 ; it is calculated from (2.3a) by 
varying the mass and keeping fixed Ao- 

In fact, the mass integration in (3.7a) can be carried out. Rewriting the right-hand side of (3.5b) as 

dU\ 1 „ dIL 

(3.8a) 



MrK + 1% T, m)) — M T K + n) ( 1 + ^ j - ^H- 



= M T(m2 + im2) _^ (m2 + ftBV _L_^n^ 



and introducing 

M{m 2 ) = f°° 2 dm' 2 M{m' 2 ) = \ f ^0^8^)9^ - m 2 ) 
M T (m 2 ) = £dm' 2 N T (m' 2 ) = J 7^T%°)% 2 ~ m ^^Zr x ^ 



we have 



M T (m 2 +U(q,T,m)) = 



M T (m 2 + 5m 2 ) + -^n 2 (T, m) 
2Ao 

so that with the boundary condition that the pressure vanishes for infinite mass 

P(T) - P(0) = M T (m 2 + 5m 2 ) - M(m 2 ) + -J- (II 2 (T, m) - U 2 (0, m)) 

2Ao 

= M T {m 2 +5m 2 ) +M{m 2 + 5m 2 ) - M{m 2 ) + \5m 2 {M T {m 2 + 5m 2 ) + M(m 2 )) 

(3-9) 

A formula similar to (3.9) has been derived previously by Amelino-Camelia and Pi [5] using the CJT 
formalism [12], though their formula does not satisfy the physically-important constraint that the 
pressure vanishes when the mass is infinite.* 

In order to highlight the interplay between thermal and quantum contributions, let us also give the 
following alternative version of the result (3.9) 

P{T) - P(0) = M T {m 2 + 5m 2 ) + ± l - ' - — M^ n ~ 1 \m 2 )(5m 2 ) n (3.10) 

n=S 

This expression makes manifest the UV finiteness** of our result, and it exhibits three different kinds 
of contribution: Mti the classical expression for the pressure of a bosonic gas of particles with mass 
squared m 2 + 5m 2 ; (5m 2 ) 2 / A, which is 0(A) in perturbation theory, essentially a thermal interaction 
contribution; and the rest, which starts at three loop order, coming from the thermal mass shift in 
zero-temperature integrals. 

* In reference [12], the effective potential is calculated, from which the pressure of our model follows by 
restricting to vanishing field expectation value and positive bare mass squared. We intend considering 
the symmetry-breaking sector of the theory in a future paper. 

** Its IR finiteness is however better seen from the original version (3.9). 
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4 Three dimensions 



When n = 3, (2.2c) and (2.9) become 



M(m ) = -- 



iV T (m 2 ) = l og (l - e - m/T ) (4.1) 

47T 

Hence the integral equation (2.10) for 5m 2 reduces to 

^-5m 2 = m- \jm 2 + 5m 2 - 2Tlog (l - exp(-y/m 2 + 5m 2 /T)) (4.2) 

A V / 

This must be solved numerically and the result plugged into (3.9). (We recall that, in 3 dimensions, 
Aq has the dimension of mass.) 



When m = 0, 



M> + £ = _21og(l-e X p(-0) (4.3) 



where 5m 2 = £ 2 T 2 . The solution for £ becomes small at high temperature, when 8ttT 3> Ao, but it 
goes to quite slowly: it is about 0.4 for 8vrT/A = 10 and 0.07 for 8vrT/A = 1000. 

In terms of £, the formula for the pressure can be written as 

P(T) - P(0) = ^ |Li 3 (e-«) + £Li 2 (e-S) - £ log(l - e"«) + (4.4) 

where Li3 and Li 2 are the tri- and di-logarithmic functions, respectively [13]. In the high-temperature 
limit, where £ — > 0, this approaches the free-field value £(3)T 3 /(27r). 

Notice that a perturbation expansion of the pressure in the m = case would be fraught with infrared 
problems. These would be encountered already at order Aq, in the two-loop graph. The formula (3.9) 
has eliminated them through resummation. The right-hand side of (3.9) is finite when m —> 0, yielding 
(4.4), but it does not have a power-series expansion in powers of 5m 2 , because the derivative of Mt 
diverges at the origin. This infrared sensitivity leads to the interesting result that at low temperature 
or strong coupling the interaction pressure becomes a constant multiple of the free-field pressure, 
independent of the coupling. For 8irT <C Ao the first term in (4.3) becomes negligible compared with 
the other two, and the solution for £ approaches 

£ = Sm/T log for T/A -»• (4.5) 

that is 5m w 0.96 T. 

Remarkably, for this value of £ the polylogarithms in (4.4) can be evaluated (see equations (1.20) and 
(6.13) of reference [13]), yielding 

P(T) - P(0) - = \[P{T) - P(0)] free for T/A - (4.6) 

07T 



We have not found any physical explanation of this surprisingly simple result nor have we been able 
to derive it without recourse to the peculiar properties of polylogarithms. 
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5 Four dimensions 

The massive case 

When in 4 we must use the integral equation (2.12) for 5m 2 , which is written in terms of the 
renormalised coupling A. With x = 5m 2 /m 2 and z = T/m, it reads 

X=[F(x,z)]- 1 

F(x, z) = ^{ (1 + x) log(l + x) - x + 8(1 + x) J™ ^ J^}_^ (5.1«) 

The function \F{x,z)\~ x is plotted against x in figure 3, for various values of T/m. Because of the 
factor 1/x in F, all the curves go to at small x. Because of the logarithm term in F, which arises 
from the (frequently neglected) contribution M in (2.11b), they go to again at large x. So for 
each, there is a maximum choice of A beyond which the equation has no real solution; for example, 
for T/m = 3 the critical value of A is about 200. For values of A below the critical value there are 
two solutions, as was found previously by Bardeen and Moshe [14]. However, only the smaller one 
is relevant, since by definition 5m 2 — ► as T — > 0; the larger solution corresponds rather to having 
a second, nonperturbative solution for the renormalized mass m as a function of the bare mass mo- 
For reasonable coupling the larger solution is in fact exponentially huge so that it would not matter 
anyway. For overcritical A, the two solutions become complex and have to be dismissed because we 
have started with the assumption of a real self-energy IT in (3.5b). 

10000 
1000 
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Figure 3: the function F 1 in (5.1a) plotted against x for various values of T/m 

Alternatively, for a given choice of A, there is a critical value of T/m beyond which there is no solution. 
For A < 10 this critical value is large and we may find it from (5.1a) approximately. The stationary 
value of the right-hand side of (5.1a) occurs when 
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_ Vk 2 ^! 
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10 




Figure 4: the function A = F 1 for m = plotted against 5m 2 /T 2 ; the heavy line marks 
where 5m would become equal to the (modulus of the) tachyon mass. 

The integral is approximately 0.8 when y/x/z = 1 and it decreases rapidly as \fxjz = 1 increases. 
Hence the stationary value is at x just greater than z 2 , and so the critical value is given by 



mexp 



That is, the critical temperature is of the order of the tachyon mass, which for A = 1 is about 10 68 m. 

The massless case 

As we have seen in section 2, the massless case requires that A — > as m — ► 0. This does not mean 
that we are driven to a trivial theory, but only that the definition of the coupling constant through 
the scattering at zero energy becomes inappropriate. Instead we have to introduce a renormalization 
scale p, and a finite coupling constant A as given by (2.3b) or (2.13). The integral equation to be 
solved for 5m 2 is now (2.14). In place of (5.1a) we have A = [F(5m 2 , /u)] _1 and the situation is similar 
to the one of the massive case with Jx replacing m. If for example we choose jl such that A(/x) = 1, the 
critical temperature above which there are no solutions is of the order of /xexp(167r 2 ). 

In the remainder of this paper, we shall compare the exact results that we have found against a 
perturbative evaluation. At a given nonzero temperature it will turn out to be useful to define the 
coupling at p, = T through 

A(/2) 



A(T) 



(5.3) 



l + \(fi)[L(P)-L(T)] 

With this choice of the renormalization scheme the solutions to A = [F(5m 2 ,fl = T)] _1 are plotted 
in figure 4. The critical value of A(T) beyond which no solutions exist turns out to be approximately 
325.5. Below this, there are two solutions for the thermal mass, but the higher solution is found in 
the region close to the tachyon scale which is marked by the heavy line in figure 4. For A(T) <C 10 2 , 
the latter is exponentially far away and this loosely defines the range of coupling where we can accept 
the theory as an effective one. 

By the way, had we chosen ji/T > 1.62 ... we would have found a seemingly different picture: for those 
values of p the critical value X(p) has moved past infinity so that there always exist two solutions for 
positive coupling; but this is only because this change of the renormalization scheme maps sufficiently 
large values of X(p = T) onto negative X(p). However, all this occurs in the region close to the tachyon 
mass scale, which we can ignore by avoiding too-large couplings. 
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We now compare the exact results one gets by a numerical evaluation of the above formulae against a 
perturbative one. This is particularly interesting in the case of a massless theory, for there ordinary 
perturbation theory runs into infrared singularities that need to be cured by resummation of the 
thermal mass. 

In perturbation theory 5m 2 /T 2 ~ A, but a naive expansion of the functions N T and Mt as a power 
series in their argument is bound to fail. However we may write the function AT T (<5m 2 ) that appears 
in the integral equation (2.14) for 5m 2 as 

T 2 ^2, 5m 

k=i 

where K\ is the modified Bessel function of the second kind. Using a Mellin transform [15] one finds 



E° 5m ( k 5m 
kf Kl 



k=i 



7r 2 7T 5m 1 / 5m Y 
~6 2T~ ~ 4 \~T~J 

2 oo 



log 



5m 



'7 



K^)s ( - ir ™ <( ** +1) (^) 



2n 



(5.4) 



With (5.4) we can solve the equation for the thermal mass perturbatively to any desired accuracy in 
A and insert the result into the pressure (3.9), which in the massless case is given by 



P(T) - P(0) = M T {5m 2 ) + \5m 2 N T (5m 2 ) + 



1 



1287T 2 



5m 4 



(5.5) 



For this, we need also 



4vr 2 , /r 2N 2vr 4 vr 2 / 5m\ 2 it ( 5m\ 3 1 f 5m x 4 
T<*r T {Sm*) =— - - - +- - +- - 



5m 3 



(^) 4 E(-> 

n=l 



1 / (5m 
4 



(^ ! c(2n | i) («HL) 



2n 



(5.6) 



The infrared divergences that appear in conventional perturbation theory without resummation of 
thermal masses have their origin in the fact that the series (5.4) and (5.6) are nonanalytic in 5m 2 at 
5m = 0. 

The first few terms of the series expansion in A are 



5m 2 



A 



A 3 / 2 



T 2 24 16ttv/6 



+ (3 - 7 - log 



p . A 2 



47rT y 3847r 2 



-(l-2 7 -21og 



A 5 / 2 



47TT' 1024^^2/3 



+ 0(A 3 ) (5.7) 



P(T) - P(0) 7T 2 



rp4 



A A 3 / 2 
+ 



90 1152 576vr^/6 



/"/ A 2 n \5/2 

" (6 " 7 " log 7^)7777^ + (3 - 2 7 - 21og -^-)- 



(6-7-lo g7 ^) 2 -30 + 



47rT y 184327r 2 
C(3)\ A 3 



36 / 294912vr 4 



+ 0(A 7 / 2 ) 



4vrT y l22887r 3 \/6 



(5.8) 
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It is simple to verify that (5.8) is independent of fx: differentiate it with respect to log/u and use 
d\/d\ogfl = A 2 /16vr 2 from (2.13). 

In the case of the N = 1 theory, the pressure has been calculated up to and including order A 5 / 2 
using hard-thermal-loop resummed perturbation theory [16]. Up to and including order A 3 / 2 , there is 
no difference between the subset of foam diagrams and the complete set of diagrams, and the results 
indeed agree. Beyond this order, there are differences, in particular there are no terms involving log A 
in the foam-diagram subset, which do occur in the full set starting at order A 5 / 2 . They come from 
the logarithmic terms in the expansion (5.4) and (5.6), which in the case of foam diagrams happen 
to combine with the log in (2.15) such that the thermal mass drops out from the arguments of the 
logarithms. 

Despite this simplification, a comparison of the above perturbative result with the exact one might 
give some hints about the convergence properties of thermal perturbation series in general. It turns 
out that these depend strongly on the ratio of jl/T. 

In figures 5-7 we juxtapose the exact and the perturbative results for the thermal mass Sm/T and the 
ratio of the pressure [P(T) — P(0)} to its ideal-gas value 7r 2 T 4 /90, including in (5.7) and (5.8) up to 10 
terms beyond the leading one. We choose various values of the renormalization scale fl, but for ease of 
comparison in each case we plot against A evaluated for fl = T through the relation (5.3). (The actual 
expansion parameter X(fl) is larger (smaller) when fl is larger (smaller) than T.) The resulting 11 
approximants are put on top of each other in order to give a visual impression of the rate (or failure) 
of convergence of the perturbative expansions in A; the exact results are indicated by dashed lines. 

When fl is very different from T, the convergence of the series deteriorates markedly. In figures 5 and 
6 the results for the thermal mass and the pressure are seen to become oscillatory for larger coupling 
when fl = WOT (the vertical lines in figures 5a and 6a are part of these oscillating results), whereas 
with fl = y^jT the perturbative results fail to improve with higher orders at roughly the same place, 
although in a less violent manner. 

With the choice fl = 2ttT, which has been advocated in reference [17] on the grounds that this is the 
mass of the first nonzero Matsubara mode, the behaviour of the series expansion can be significantly 
improved, although for the first few approximants the results are slightly worse than at fx = T, see 
figure 7. 

Attempting to improve perturbation theory by putting to zero one of the A 2 -terms in (5.7) or (5.8) 
gives much larger fi's and rather bad convergence properties. The optimal choice of fl seems to be 
around fl = 47rexp(— 7)T, which absorbs all 7's and log(47r)'s. This is in fact rather close to the 
choice fl = 2ttT of reference [17]. Notice that the origin of the 7's in (5.7) and (5.8) is entirely from 
the high-temperature expansions (5.4) and (5.6); those appearing in dimensional regularization have 
already been absorbed in fl. However, without the restriction to foam diagrams, the 7's would not have 
the same coefficients as the logs, so it is not clear whether this result for the optimal renormalization 
point could be a general one. But it does confirm the expectation [17] that the optimal renormalization 
scheme is to be found around 2irT rather than T. 

Another noteworthy observation is that the rate of convergence is markedly slower for the perturbation 
series of the pressure than it is for the thermal mass. Since this loss of accuracy comes from having 
inserted a perturbative result for Sm/T into the series (5.6) and truncated at a given order in the 
coupling, a certain improvement would be simply to refrain from doing a high-temperature expansion 
of the integrals that appear in the expression for the pressure. The quality of the perturbation series for 
the pressure is then the same as that of the thermal mass. It would be interesting to see whether this 
could be implemented in QCD to ameliorate the frustratingly bad apparent convergence of resummed 
perturbation theory for the QCD pressure, which has been calculated up to order g 5 recently [18]. 
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(a) 




(b) 



A*(T) 




Figure 5: A comparison of the perturbative results for Sm/T as a function of A X / 2 (T) up 
to 10th order for different choices of the renormalization scale: a) p = 100T, b) p = T, c) 
p, = y^qT. In a) the "vertical" lines are part of two of the curves on the left. 
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Figure 7: The perturbative results for Sm/T and P = [P(T) — P(0)]/ 7r 90 for p = 2ttT. 
Compared to figures 5b and 6b, the rate of convergence becomes much more rapid after the 
first few approximations, for which p = T is slightly favoured. 

A rather simple reorganization of perturbation theory is brought about by replacing truncated power 
series expansions by perturbatively equivalent Pade approximants [19], which amounts to replacing 



F n = c + cig 1 + 



c n g 



\p,n-p] 



Co + Qiff 1 + ■ ■ ■ + a p g p 
1 + cip+ig 1 + . . . + a n g 



n—p 



+ 0(g n+1 ) 



(5.9) 



In reference [20], this possibility has been studied in the context of thermal perturbation theory by 
testing for an unphysical dependence on the renormalization scale. Our results allow us to investigate 
the quality of Pade approximants by comparing directly with truncated power series and the exact 
result. 

In figure 8 we have replaced the various nth-order results that gave rise to the curves displayed in 
figure 6b by [n/2,n/2] Pade approximants up to n = 8 (for odd n we rounded off the second entry 
at the expense of the first) and using p, = T. It turns out that we find a spectacular improvement of 
convergence up to really high values of the coupling which in our theory is bounded by the requirement 
to be sufficiently below the critical value, A X / 2 (T) <C 18. The lower line in figure 8 is the Pade 
approximant [0, 2] ([1, 1] does not exist), which is only marginally better than its perturbative analogue 
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Figure 8: As in figure 6b, but with the perturbative power series replaced by Pade approxi- 
mants [0, 2], [1, 2], [2, 2], [2, 3], ... , [4, 4]. They quickly converge to the exact result with the 
exception of [3, 3], which has a pole beyond which this approximant seems off by a constant. 

[2,0], but already [1,2] (given by the next line upwards) is a very good approximation. The higher 
approximants approach the exact result from above and are extremely accurate with the exception 
of [3, 3] which has a pole at A 1 / 2 « 9.5. For A 1 / 2 sufficiently smaller than that, this approximant is 
in fact quite good, but after the pole has been encountered it seems to be off by a constant* and its 
quality is inferior to most of the lower-order approximants. 

Curiously enough, with the choice p = 2itT which has led to more rapid convergence in the case of 
truncated power series, the Pade approximants are not as good as in figure 8, although they are still 
a great improvement for A 1 / 2 < 8. What happens is that all of the approximants run into poles in the 
range of A considered. The same holds true for larger p, whereas for p <C T the quality of the Pade 
approximations decreases, too, but not as rapidly. 

All in all it appears that Pade approximants can give vast improvements of a truncated power series 
expansion unless the rational functions used as approximants develop poles. While there is no real 
theoretical explanation for the superiority of Pade approximants, it seems that their main advantage 
is that (in the absence of poles) the latter do not blow up at larger coupling as quickly as the corre- 
sponding truncated power series inevitably do. Since the exact result behaves rather unspectacularly, 
the odds seem to be in favour of Pade approximants. 

Unfortunately, in QCD Pade approximants turn out to lead to less impressive improvements [20]. In 
figure 9a the perturbative result for QCD with 3 flavours is given, which shows that (resummed) 
perturbation theory is useful only up to g{T) « 1/2. But a real quark-gluon plasma as one hopes 
to produce in heavy-ion collisions has rather g(T) ~ 2, where the perturbative results are completely 
inconclusive. The corresponding Pade approximants are rendered in figure 9b. They seem to give 
some improvement of convergence, extending the allowed range of coupling to perhaps g(T) f» 1, but 
this appears to break down before reaching g(T) = 2. 

So ultimately one would have to find a different expansion scheme that does not involve truncated 
series in the coupling if one wants to cover more-strongly-coupled theories. Recently an interesting 
attempt towards an alternative perturbative scheme has been made in the example of a scalar theory 
in reference [21] using the numerical solution of an approximation to the gap equation in a loop 

* In reference [20] poles in Pade approximants were simply subtracted. Our observations imply that 
this works nicely only when the coupling is such that one is below the point where a pole arises. 
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Figure 9: a) The perturbative results for the pressure of QCD(Nf = 3) up to order g 5 (full 
line). Short, medium, and long dashes give the results up to g 2 , g 3 , and g 4 , respectively, b) 
The corresponding Pade approximants [0,2], [1,2], [2,2], and [3,2]. 



expansion of the pressure, which however required an ad hoc treatment of uncancelled ultraviolet 
divergences. 

In this context we observe that while the perturbative results are satisfactory only for small coupling, 
the high-temperature series (5.4) and (5.6) have excellent convergence properties, with convergence 
radius 5m /T = 2tt. In our model, the exact result for the thermal mass remains sufficiently smaller 
than 2ttT even for extremely large coupling A. Only when we re-expanded in A did the convergence 
become bad. In order to highlight this, figure 10a shows the solutions to the thermal mass equation 
when it is truncated at the leading or next-to- leading term in the expansion in 5m /T, rather than 
being expanded in A. Including the next power of 5m 2 /T 2 , together with the logarithm, turns out 
to produce a result which is virtually indistinguishable from the exact one for the whole range of A 
bounded by the requirement that the tachyon mass remains larger than anything else. The same holds 
true for the pressure (figure 10b), with even smaller deviations from the exact result. 

In the particularly interesting case of gauge theories it is crucial to have a consistent expansion scheme 
with a well-defined expansion parameter in order to be able to retain gauge fixing independence. If 
in QCD it would be possible to reorganize perturbation theory as a series in 5m /T rather than g, 
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Figure 10: a) Solutions of the thermal mass equation when the latter is truncated at leading 
order in 5m /T (top curve) and next-to-leading order (bottom curve). All higher approxima- 
tions are virtually indistinguishable from the exact solution, and they collectively form the 
middle line, which broadens only at the highest values of A. b) Analogously for the pressure. 

figure 10 suggests that this should lead to dramatic improvements. 
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